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ABSTRACT 



We report high-resolution (0.07 arcsec) near-infrared polarized intensity im- 
ages of the circumstellar disk around the star 2MASS J16042165-2130284 ob- 
tained with HiCIAO mounted on the Subaru 8.2 m telescope. We present our 
if -band data, which clearly exhibits a resolved, face-on disk with a large inner 
hole for the first time at infrared wavelengths. We detect the centrosymmetric 
polarization pattern in the circumstellar material as has been observed in other 
disks. Elliptical fitting gives the semimajor axis, semiminor axis, and position 
angle (P.A.) of the disk as 63 AU, 62 AU, and -14 °, respectively. The disk is 
asymmetric, with one dip located at P.A.s of ~85°. Our observed disk size agrees 
well with a previous study of dust and CO emission at submillimeter wavelength 
with Submillimeter Array. Hence, the near-infrared light is interpreted as scat- 
tered light reflected from the inner edge of the disk. Our observations also detect 
an elongated arc (50 AU) extending over the disk inner hole. It emanates at the 
inner edge of the western side of the disk, extending inward first, then curving to 
the northeast. We discuss the possibility that the inner hole, the dip, and the arc 
that we have observed may be related to the existence of unseen bodies within 
the disk. 

Subject headings: stars: pre-main sequence — planetary systems — protoplanetary 
disks — techniques: polarimetric 



-5 - 



Introduction 



Planets are believed to form in protoplanetary disks. To study the structures of 
protoplanetary disks, multi-wavelength observations have targeted such disks around young 



Scor 3 


ius regions (e.g.. 


Kenvon & Hartmann 


1995; 


Preibisch & Zinnecker 


1999; 


Fur lan et al. 


2006: 


Andrews et al. 


2009; 


Dahm & Carpenter 


2009; 


McClure 


2009; 


Andrews et al. 


2011 




Mathews et alJboia) 


Among such nearby star forming regions, Upper Scorpius, which 





is located at a distance of 145 pc f ide Zeeuw et al. 



19991 ). provides a 



particularly in t erms of its age. While the age is only ~5-10 Myr fjPecaut et al 



Song et al. 



unique environme nt 



2012 



20121 ). the lack of dense molecular material and embedded young stellar objects 
indicate that the star formation process in Upper Scorpius is over. The area is essentially 
free of dense gas and dust clouds, and the association members show only moderate 
extinctions (Ay < 2 mag). In fact, there are only two K and M stars (1.5%), including the 
star discussed here, in the Upper Scorpius region which ex hibit a if-band ex cess, although 
about half of the stars in Taurus exhibit a f^-band excess ( jStrom et al.lll989l ). This implies 
that there are very few stars still hosting protoplanetary disks in Upper Scorpius star 
forming region. 

Here, we introduce our target 2 MASS J16042165-2130284 (hereafter J1604-2130) in 
Upper Scorpius. J1604-2130 was first identified as a member of th e Upper Scorpius OB 



Preibisch fc Zinnecker 



association in the spectroscopic survey of X-ray selected sources of 
( 119991 ). who found: stellar age 3.7 Myr, spectral type K2, stellar mass 1 M , 
logT c fj[K]=3.658, and \og[L/L ©] = —0.118. The s pectral energy distributio n (SED) of 



J1604-2130 has the features of a transitional disk ( iDahm fc Carpenter 



20091 ). in which 



1 Based on data collected at the Subaru Telescope, which is operated by the National 
Astronomical Observatory of Japan. 
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small dust grains have been partly cleared 



larger radii. Recently, 



Mathews et al 



rom the inner disk but are still abundant at 



(120121 ) presented the Sub-Millimeter Array (SMA) 



880 /im images, which directly resolved a face-on protoplanetary disk around J1604-2130, 
and derived a dust mass of 0.1 Mj up . The availability of radio data was one advantage to 
selecting this target, however other advantages abound and are as follows: 1. J1604-2130 
is located in a nearby star forming region, Upper Sco; 2. it hosts one of the largest inner 
holes in transitional disks discovered to date; 3. it hosts t he most massive disk in Upper 



Scorpius; 4. the inclination 6° ± 1.5° f lMathews et al. 



20121 ) is very close to face-on. 



Infrared photometry and spectroscopy with Spitzer has shown that less than ~10% of 



stars with disks s 



(e.g., 



row some degree of inner disk depletion yet still r etain massive outer disks 



Cieza et al 



2010 



Muzerolle et al. 



2010 



Luhman et al 



20111 ). Although there have 



been several transitional disks detected with radio observations with S MA around single 



young solar mass sources (e.g. 



Pietu et al. 



2006 



Andrews et al. 



of transitional disks have been detected in th e near-infrared (e.g 



Thalmann et al. 



2010 



Hashimoto et al. 



2011). only a handful 



Fukagawa et al. 



2006 



20HI). This is mainly because it is difficult to 



resolve structures with only few 10 AU scale close to a bright central star. Furthermore, 
structures inside the hole of disks have rarely been directly imaged or resolved to date. 
To gain a better understanding of the morphological structures at radii where planets are 
expected to form, we conducted high-resolution near-infrared polarimetric observations of 
J1604-2130. 



2. Observations and Data Reduction 



We carried out polarimetry in ff-b and (1.6 /xm) toward J1604-2130 using t he high 
resolution imaging instrument HiCIAO (ITamura et al.l 120061 ; iHodapp et al.l |2006[ ) with a 
dual-beam polarimeter at the Subaru 8.2 m Telescope on 2012 April 11. These observations 



are part of the ongoing high-contrast imaging survey (SEEDS; iTamural 120091 ). The 
polarimetric observation mode acquires o-rays and e-rays simultaneously, and images a 
field of view of 10"x20" with a pixel scale of 9.53 mas/pixel. J1604-2130 was observed 
without an occulting mask in order to image the inner most region around the central star. 
The exposures were performed at four position angles (P.A.s) of the half-wave plate, with 
a sequence of P.A. = 0°, 45°, 22.5°, and 67.5° to measure the Stokes parameters. The 
integration time per wave plate position was 15 sec and we obtai ned twenty-five waveplate 



cycles. The adaptive optics system (AO 188; 



Hayano et al. 



20101 ) provides a diffraction 



limited and mostly stable stellar PSF. The total integration time of the polarization 
intensity (hereafter PI ) image was 750 sec after removing low quality images with large 
FWHMs by careful inspections of the stellar PSF. 

The Image Reduction and Analysis Facility (IRAFg) software wa s used for all data 



reduct ion. The polarimetric d ata reduction technique is described by 



( 1201 lh and 



Muto et al 



Hashimoto et al. 



(120121 ) . using the standard approach for polarimetric differential 
imaging (Hinkley et al. 2009). By subtracting two images of extraordinary and ordinary 
rays at each wave plate position, we obtained +Q, —Q, +U, and —U images, from which 
2Q and 2U images were made by another subtraction to eliminate remaining aberration. 
The PI was then calculated by PI = \/ Q 2 + U 2 . The i nstrumenta 



at the Nasmyth platform was corrected for by following 



Joos et al 



polarization of HiCIAO 



((20081). 



2 IRAF is distributed by National Optical Astronomy Observatory, which is operated by 
the Association of Universities for Research in Astronomy, Inc., under cooperative agreement 
with the National Science Foundation. 
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3. Results 

The if-band PI image of J1604-2130 is presented in Figure 1(a). The emission 
arises from dust particles mixed with gas in the circumstellar structures, scattering the 
stellar light. Details of the near-infrared disk around J1604-2130 have for the first time 
emerged at 0.07 arcsec. The disk has a dip D, located around a position angle ~85°, 
making this transitional disk asymmetric. Inside the bright disk, we see a region where 
the surface brightness drops by a factor of ~5 compared with the disk. We call this inner 
lower-brightness region the "hole" throughout the remainder of this paper. Figure 1(b) 
shows if-band polarization vectors. The disk with its centrosymmetric polarization pattern 



surrounds the central source J1604-2130 as previous 
technique of near- infrared imaging polarimetry (e.g., 



y observed in ot 



rer di sks using the 



Hashimoto et al. 



201lh . 



Elliptical fitting was performed to measure several disk parameters including the 
angular separations of the major and minor axes, the position angle, inclination, and the 
position of the ellipse center. We first measured disk radial profiles in 5° position angle 
increments, and then extracted coordinates of the brightest peak area. Those coordinates 
were used to fit an ellipse. The nonlinear least-squares Levenburg-Marquardt algorithm 
with the above five free parameters was performed and its result is displayed in Table 1. 
Its shape is superimposed on Figure 1(a). The derived semimajor and semiminor axes are 
called "NIR peak radii" in this paper. The geometric center of the disk is consistent with 
the position of the star within errors. The disk position angle has a relatively large error 
because the disk is close to circular. 

Figure [2] shows the azimuthal surface brightness at r = 33, 63, 145 AU, which are 
illustrated in figure 1(c), with the position angle measured from north to east. The error 
bars shown in figure [2] represent the calculated standard deviation. The left panel of figure 
[3] shows the radial surface brightness profile of J1604-2130 along the minor axis. The error 
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bars shown in figure [3] represent the calculated standard deviation. In the north, the surface 
brightness along the minor axis increases as r 3 - 6±0 2 from 40 to 62 AU, and decreases as 
r -4.7±o.i f rom 74 ^ 124 AU. In the south, the surface brightness along the minor axis 
increases as r 3 ' 0±0,4 from 40 to 62 AU, and decreases as r~ 4 - 7=t0,2 from 74 to 124 AU. The 
right panel of Figure [3] shows the radial profile of the surface brightness along the major 
axis. In the east, the surface brightness along the major axis increases as r 2,6±0 - 2 from 40 
to 62 AU, and decreases as r~ 4 - 7±ai from 74 to 124 AU. In the west, the surface brightness 
along the major axis increases as r L9±a2 from 40 to 62 AU, and decreases as r _4,0±0 - 2 from 
74 to 124 AU. The profiles along both the major and minor axes show a change of slope 
beyond 63 AU, thus reconfirming the NIR peak radius to be 63 AU. Along all four axes, the 
counts drop to half of the peak counts at around 50 AU from the central star. The width of 
the brightest region where counts exceeds half of the peak counts in the disk is around 30 
AU. 

A salient feature inside the hole is an elongated arc extending over the disk inner hole. 
It begins at the inner edge of the western side of the disk, extending inward first, then 
curving to the northeast. The whole arc extends at least 50 AU. This arc is equivalent in 
brightness to the inner edge of the disk and to the dip D region. This arc consists of a 
real feature and is not an artifact because in contrast to the inner hole errors (shown as 
squares in Figure |2]), the surface brightness difference between the arc and the other regions 
in the inner hole has a 6a significance, whereas the other directions show small fluctuations 
in emission. Furthermore, the polarization vectors in the region of the arc mostly face 
the central star as well as the region of the bright disk, indicating that the arc is also 
illuminated by central star and is physically connected to inner edge of the disk. 

In the same way, the dip D is a real dip feature and is not an artifact because the 
difference in surface brightness between D and the rest of the disk is statistically significant 



- 10 - 



and not within errorbars. In addition, the polarization vectors in the region of the dip 
mostly face the central star. 

There appears to be another dip located around a position angle ~255°, which is not 
azimuthally symmetric with D in terms of location and surface brightness. However, this 
another dip is not significant outside of the uncertainties of the observations, thus follow up 
future observations are required to confirm its presence. 



4. Discussion 



4.1. Comparison with the SMA Imaging 



Mathews et al.l (120121 ) present SMA 880 /im continuum and CO (3-2) line observations 
of this source with a resolution of 0".34 and resolve the disk (Figure 1(d)). They determine 
the surface density distribution of the disk with inner cavity whose radius is 72 AU based 
on the model fitting to the visibility at 880 /im. They also suggest that sparse matter 
still remains inside the cavity (20 - 72 AU) because of the dust emission indicated by 



mid- infrared (ICarpenter et al. 



20091 ) excesses and their model fitting. The NIR peak disk 
radius of 63 AU derived from our if -band image roughly agrees or may be slightly smaller 
than their measurement and can be explained by scattered light reflected at the inner edge 
of the disk. 



Based on their CO peak velocity map, Mathews et al.l ( 120121 ) derive the position angles 
of the major and minor axes to be 85° ± 10° and —5° ± 10°, respectively. Their derive 
disk parameters, including the position angle and inclination, are in a good agreement with 
values derived from our near infrared observations. 



Mathews et al. 



(120121 ) point out that the relatively high mass of the gas-dust disk, its 



sharp inner edge, and the lack of gas accretion onto the star together point to a dynamical 
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origin for the inne r hole. No compani ons have been found down to brown dwarf masses 



(Kraus et al. 


2008; 


Ireland et al. 


2011) 



origin of the hole can be attributed to the unseen presence of one or several giant planets 
within about a 40 AU. 



The width and depth of the hole depend on the 



disk thickness (e.g., 



Takeuchi et al 



1996 



Crida et al. 



plane t mass, disk viscosity, and 



20061 ). As discussed by Dong et 



al. (2012)b, the radius of the inner hole is on the order of ~ 60 - 70 AU. There have been 
extensive studies on the process of gap opening by a planet em bedded in a disk (e.g., 



Lin fc Papaloizou 



1986 



Takeuchi et al 



1996 



Crida et al. 



20061 ). A rough estimate of the 



upper limit of the size of the gap that one planet can make may be given by the location of 
the lowest order Lindblad resonances, which, in the case of the outer resona nce (the outer 



edge of the gap), is at ~1.6 times the orbital radius of the embedded planet (jTakeuchi et al. 



19961 ). If this is the case, the existence of the inner hole might indicate the existe nce of at 



Crida et al. 



20061 ). and the 



least one planet at ~40 - 50 AU from the central star (see also 
inner edge of the hole would be expected at ~ 25 — 30 AU, which is inside of our saturation 
radius and is therefore unresolved in this study. 

If multiple planets reside within the hole, the hol e width can be w ider and the inner 



201 ll ). To understand the 



disk may be much smaller than the single planet case (iZhu et al. 
properties of the whole of this object, it is important to give constraints on the inner disk 
structure using detailed SED modeling (see Dong et al. (2012)b) and/or high-resolution 
observations by ALMA. Our observations do not show a clear inner disk edge, suggesting 
that the inner disk radius must be inside of the saturation radius of our observations. 
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4.2. Non-axisymmetric Features 



We list a few possibilities for origin of the two non-axisymmetric features, the dip and 
th e arc since it is difficu lt to draw firm conclusions at this stage. According to Figure 6.(c) 



in 



Mathews et al 



( 120121 ) . the only structures that are not reproduced by the axisymmetric 
disk model within 2a error are the south-east component; the north dip is consistent with 
the axisymmetric disk model convolved with the elliptical beam. We, ther efore, do not 



elabor ate on the north dip on the sub-mm disk in the following discussions. iMathews et al. 



(120121 ) observe that the millimeter continuum emission has a peak at 0".5 southeast from 
the star in the 880 /im image. A similar distribution is also seen in the total integrated 
intensity distribution of CO (3-2) emission. Our if-band image and Figure [3] show that the 
southeast of the disk is brighter than the east and the west, despite the fact that we observe 
the scattered light originating at the disk surface. This may be explained in a way that the 
more materials exist, the thicker the disk is, resu lting in a higher grazing angle of the stellar 



light at the scattering surface (e.g. 



Muto 



201lh . In our near- infrared image, the north of 



the disk is also bright. The brightness of the disk in the north, together with sub-mm 
observations suggesting material accumulation in the south, infers that the northern portion 
of the disk is the near side. This assumes that forward scattering dominates. If this 
picture of the disk orientation is correct, the disk rotation is counter-clockwise direction. 
This geometry might also explain the origins of the non-axisymmetric structures at D. D 
reside at the interface between where the material accumulated (south) and where forward 
scattering is significant (north). However, it should be noted that since the disk is nearly 
face-on, significant forward scattering excess and precise prescriptions of the phase function 
should be required, and our discussion crucially depends on the details of the optical 
properties of the dust particles. It should also be n oted that the materia l concentration in 



the south is within 2a of the current sub-mm data (IMathews et al. 



20121 ) . The detection of 



the non-axisymmetric structures at sub-mm wavelengths is needed to examine the details 
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of the disk geometry. 



Another possibility for the origin of the dip D is that there is a massive gravitating 
object at the disk's midplane. It could be a planet (which is not necessarily the same as the 
one discussed in the previous section as the origin of the inner hole), a clump formed by 
hydrodynamic instability (e.g., gravitational instability), or something else. Such objects 



can pull the materials from the sur face toward the dis 



shadow in the scattered light image (I Jang- Condell 



i's midplane, thereby producing a 



20091 ). This might produce features such 



as D. Such objects have sizes smaller than the SMA beam size, therefore cannot be resolved 
with the SMA. 

For the arc structure in the west, the morphology is similar to spiral density waves in 
the disk. An arguably uniqueness is that this arc extends over the "disk inner hole" region. 
It also resembles to spiral density waves in terms of the pitch angle. Spiral waves can be 
produced by dynam ical processes in the di sk such as the turbulence or the planet-disk 



interaction (see e.g., 



Papaloizou et al. 



20071 ) . If this feature is a spiral wave, its morphology 
indicates that the disk rotation is counter-clockwise because spirals are in general trailing 
features. This rotation direction is consistent with the picture presented at the beginning 
of this section. Last, we note t hat the arc structure shows morphological similarity to the 



"planet shadows" discussed by 



Jang-Condelll (120091 ) . One possible way to distinguish the 



different scenarios is to observe the structures of the disk's midplane at longer wavelength. 
For instance, if a density enhancement is observed at the location of D, this might indicate 
the existence of a massive body there. If we see an arc structure in the sub-mm as well, 
then it is likely a density wave. 

Another possible way to reveal the nature of the arc may be the contemporaneous 
SED and scattered light observations. Previous studies indicate that substantial variability 
has also been seen in MIR SED and has been linked to variable structures including arcs 
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or rims of the c 


isk ( 


Dahm & Carpenter 


2009: 


Wisniewski et al. 


2008: 


Muzerolle et al. 


2009; 


Espaillat et al. 


2011 


Flahertv et al. 


201 


2). For J1604-2130 case, the location of the arc and 



the rim of the disk wall is generally consistent with the radial location in the disk which 
possibly contributes to the SED variability. 
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Fig. 1. — if-band HiCIAO images of J1604-2130. The saturated central area (radius = 
0."2) is masked in black, (a): The PI image of J1604-2130. The field of view (FOV) is 
2". 9 x 2". 9. The unit of the color bar is mJy/arcsec 2 . The light blue ellipse and plus sign are 
the best fit result of our elliptical disk model and the ellipse center, (b): if-band polarization 
vectors superposed on the PI image. The vector directions indicate angles of polarization. 
The plotted vectors are based on 7 [pixel] x 7 [pixel] binning corresponding to the spatial 
resolution. The FOV is 2".0 x 2".0. The vector's lengths are arbitrary. (c):Red(r= 145 AU), 
brown(r = 63 AU), and yellow(r = 33 AU) circle s, corresponding figu re El superimposed on 
the PI image. (d):SMA 880 /im continuum map (IMathews et al.ll2012l ) superimposed on the 
PI image. White color contours indicate 2, 3, 6, 9, and 12<r intensity (la = 1.3 mJy/beam). 
The ~ 0."5 x 0."3 beam of SMA is shown in the bottom right. 
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Fig. 2. — Azimuthal surface brightness profile at r = 33 AU (inner hole region), r = 63 AU 
(inner disk region), and r— 145 AU (outer region) with position angle measured from north 
to east. The radius 33 AU is measured as it is the brightest part of the arc. 
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Fig. 3. — Left Panel:Radial profile along the minor axis. NW and SE radial profile are 
averaged over -25° <P.A. < 5° and 155° <P.A. < f 85°, respectively. This is because the disk 
position angle derived from our near-infrared image and peak velocity of the submm CO 
emission (—5° ± 10°) are 9° different. Therefore, we centered our 30° wide radial profiles on 
the average of our near-infrared and submm line emission. Right PanekRadial profile along 
the major axis. In the same way, NE and SW radial profiles are averaged over 65° <P.A. 
< 95° and 245° <P.A. < 275°, respectively. 
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Table 1: The result of the ellipse fitting for the disk of J1604-2130 

Semimajor axis 63.4 ± 1.0 [AU] 

Semiminor axis 62.4 ± 1.0 [AU] 

Ellipse center (u,v) a (-1.1 ± 1.1, 1.5 ± 1.1) [pixel] 

Position angle — 14 ± 11° 

Inclination 10.2 ± 10.1° 



"Central position (0, 0) corresponds to the stellar position. 



